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Abstract. Our goal is to identify the type and number of static equilibrium 
points of solids arising from fine, equidistant n-discrctrizations of smooth, con- 
vex surfaces. We assume uniform gravity and a frictionless, horizontal, pla- 
nar support. We show that as n approaches infinity these numbers fluctuate 
around specific values which we call the imaginary equilibrium indices asso- 
ciated with the approximated smooth surface. We derive simple formulae for 
these numbers in terms of the principal curvatures and the radial distances 
of the equilibrium points of the solid from its center of gravity. Our results 
are illustrated on a discretized ellipsoid and match well the observations on 
natural pebble surfaces. 



1. Introduction 

The study of equilibria of rigid bodies was initiated by Archimedes [T] ; his results 
have been used in naval design even in the 18th century (cf. 2J). Archimedes' main 
concern was the number of the stable balance points of the body. 

Static equilibria of convex bodies correspond to the singularities of the gradient 
vector field characterizing their surface. Modern, global theory of the generic sin- 
gularities of smooth manifolds appears to start with the papers of Cayley [4j and 
Maxwell [Mj who independently introduced topological ideas into this field yielding 
results on the global number of stationary points. These ideas have been further 
generalized by Poincare and Hopf, leading to the Poincare-Hopf Theorem [3J on 
topological invariants. If applied to generic, convex bodies, represented by gradient 
fields defined on the sphere, this theorem states that the number S of 'sinks' (stable 
equilibria), the number U of 'sources' (unstable equilibria) and the number N of 
saddles always satisfy the equation 

(1) S + U-N = 2. 

This formula, the so-called Poincare-Hopf formula can be regarded as a generaliza- 
tion of the well-known Euler's formula |12j for convex polyhedra. 

Static equilibria of polyhedra have also been investigated; in particular, the 
mininal number of equilibria attracted substantial interest. Monostatic polyhedra 
(i.e. polyhedra with just 5 = 1 stable equilibrium point) have been studied in |13j . 
0,0 and 0. 
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The total number T of equlibria (T = S + U + N) has also been in the focus of 
research. In planar, homogeneous, convex bodies (rolling along their circumference 
on a horizontal support), we have T > 4 [5]. However, convex homogeneous objects 
with T = 2 exist in the three-dimensional space (cf. |17|). Zamfirescu [18] showed 
that for typical convex bodies, T is infinite. While typical convex bodies are nei- 
ther smooth objects, nor are they polyhedral surfaces, Zamfirescu's result strongly 
suggests that equilibria in abundant numbers may occur in physically relevant sce- 
narios. 

This is indeed the case if we study the surfaces of natural pebbles which, while 
rolling on a horizontal plane, are supported on their convex hull Their convex 
hull is well approximated by a many-faceted polyhedron P, on which, by studying 
its detailed 3D scanned images one can observe large numbers of adjacent equilibria 
in strongly localized flocks. If we approximate the polyhedron P by a sufficiently 
smooth surface M, then we can see that the flocks of equilibria on P appear in the 
close vicinity of the (isolated) equilibrium points of M (cf. Figure [4]). 

In this paper we seek a mathematical justification for this observation. We study 
the inverse phenomenon: namely, we seek the numbers and types of static equilib- 
rium points of the families of polyhedra P n arising as equidistant n-discretizations 
on an increasingly refined [u, u]-grid of a smooth surface M with T generic equi- 
librium points, denoted by m; (i = 1,2, . . . T). As n — > oo, P n —> M and we find 
that the diameter of each of the T flocks on P n (appearing around m,) shrink and 
approach zero. However, we also find that inside a fixed, 2k x 2k rectangular grid 
domain (centered at m,), the numbers S™, U™, NJ 1 of equilibria in each flock fluctu- 
ate around specific values S* ,U* and N* that are independent of the mesh size and 
the paramctrization of the surface. We call these quantities the imaginary equilib- 
rium indices associated with m*. We may eliminate the fluctuation of S*,U*,N* 
by averaging over meshes in random positions (with uniform distributions) and in 
Theorem [2] of Section [4j we obtain the following simple formulae in terms of the 
principal curvatures K\,i,H2.i < and the distance pi of rrii from the center of 
gravity: 

(2) S* = di, U* = Ki ti K 2 ,ipjdi, N* = + K2,i)p l d i 

where di = l/\(Ki t iPi + 1)(k2,iPi + We also remark that due to (|2| we have 



and thus, by summation over all equilibrium points rrii of M , for the full polyhedron 
P n , the Poincare-Hopf formula is satisfied. 

Since imaginary equilibrium indices are defined via local quantities associated 
with mi, we may interpret them as describing properties of infinitely fine discretiza- 
tions. Nevertheless, they also provide approximate prediction for the number and 
type of localized equilibrium points for finite, dense meshes. Beyond those numbers, 
we also identify the spatial patterns associated with equilibria on fine discretizations 
and illustrate these patterns on an ellipsoid (cf. Figure [T]). We also note that these 
numbers match well with observations on scanned pebble surfaces, as we point out 
in Section [6j Remark [7] and illustrate this in Figure [4j 

In Section [2] we introduce the notations and concepts used in our investigation. 
We examine the properties of the equilibrium points of curves in R 2 in Section pi 
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Our results about those of a surface in R 3 are presented in Section [3] In Section [5] 
we describe a method which, under a slightly special condition, finds a more direct 
connection between imaginary equilibrium indices and the numbers of equilibrium 
points of the approximating surface than the one in the previous sections. Finally, 
in Section [6j we collect our remarks, and propose some open problems. 

2. Preliminaries 

In this paper, we identify points with their position vectors. For points p, q £ 
M. d , we denote the closed segment with endpoints p, q by [p,q\. Furthermore, we 
denote the standard inner product of p and q by (p, q) , and the Euclidean norm 
of p by \\p\\ = \J (p, p) . The convex hull of a set S is denoted by convS", or, if 
S = {p\,p2, ■ ■ ■ ,Pk} is finite, then possibly by \pi,P2, ■ ■ ■ ,Pk]- We denote the origin 
by o. 

Since we intend to examine equilibrium points of nonconvex and/or nonsmooth 
surfaces, we introduce the following types of equilibrium points, which we examine 
for d = 2 and d = 3. 

Definition 1. Let H C K d be a C 1 -class hypersurface (i.e. a compact, 1-codi- 
mensional C 1 -submanifold in M. d ), and let p £ R d \ H. If the (unique) tangent 
hyperplane of H at m £ H is orthogonal to m — p, we say that m is an equilibrium 
point of H relative to p. 

Definition 2. Let P c K d be a polyhedral hypersurface (i.e. a compact, 1-co- 
dimensional C°-submanifold in R d which is also the polyhedron of a polyhedral 
complex), and let p £ M. d \ P. If, m £ P has a neighborhood V such that the 
hyperplane H, orthogonal to m — p and passing through m, supports conv({p} U 
(P n V)), we say that m is an equilibrium point of P relative to p. 

For simplicity, we assume that our reference point is the origin o. We deal only 
with generic equilibrium points of C 2 -hypersurfaces (cf. |15j); that is, we assume 
that at these points the Jacobian of the Euclidean norm function is regular. In 
other words, if p denotes the distance of o and the equilibrium point m, we assume 
that k 7^ — i, where k is the (signed) curvature of the curve at m for d = 2, and 
any of the two fundamental curvatures of the surface at m if d — 3. 

3. The equilibria of plane curves 

Throughout this section, we assume that the curve under investigation satisfies 
the C 3 differentiability property, and has exactly one equilibrium point m. Note 
that as a plane curve is a one-dimensional submanifold of M 2 , there is a neighbor- 
hood of m in which the examined curve is given as a simple r : [n, t 2 ] — > M 2 three 
times continuously differentiable curve. Since for sufficiently fine discretization, all 
the equilibrium points of the approximating polygon are located in this neighbor- 
hood of m, we may assume that the curve is given in the above form. For simplicity, 
we may assume that m = r(0) = (0, p) with £ {r\,T2) and p > 0, and that the 
tangent line at r(0) is horizontal; or in other words, that y(0) = 0. Furthermore, 
we may assume that f(0) > 0. We denote the signed curvature of r at m by n. 
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Let F n denote the n-segment equidistant partition of [ti,T2]. If ^ T2 is a 

divison point of F n , then we introduce the notation pf = r ( ^ T2 ^~ Tl ^ . Let P n be 

the polygonal curve with edges of the form [p2-ijP?]- Note that then P n is the 
approximating polygonal curve defined by the equidistant partition F n , and, fur- 
thermore, the vertices of P n are labeled in a way that the indices are not necessarily 
positive integers, but real numbers that are congruent mod 1. These numbers can 
be written in the form — h k, where k = 0, 1, . . . , n. 

Let U n (K) denote the number of the equilibrium points of P n at the vertices 
p™ of P™ satisfying |i| < K. Similarly, let S n (K) be the number of the equilibrium 
points of P n lying in the relative interiors of edges [p",p" + i] satisfying |i| < K. 
Our aim is to determine the values of U n {K) and S n (K ) for 'large' values of n and 
K . We assume that the equilibrium points of P n are generic; that is, that for any 
equilibrium at a vertex p™, the vector p™ is orthogonal to neither p™ — p"_i, nor to 
Pi ~ Pi+v The main result of this section is the following. 

Theorem 1. Let K + K(n, p) be sufficiently large. Then there is a value of n e such 
that for every n > n e , we have 



U n {K) 

and 

S n {K) 





\pn\ 


or 


\pn\ 




_\ P K + 


1|_ 


_\ P K + 


1|J 




1 




or 


1 






_\pK + 


11 j 


_\PK + 


1|J 



Proof. We use the notations introduced in the preliminary part of this section, and 
for simplicity, instead of p" we write only pi . Let A„ = T2 ~ Tl , which yields that 
Pi = r(iA n ). 

First, we consider the relative interiors of the edges of P n . Our key observation 
is that [pi, pi+i] contains an equilibrium point if, and only if the following holds: 

(3) (Pt,Pi+i - Pi) < < (p t+1 ,p l+1 - pi). 

To find the equilibrium points, we use the second degree Taylor polynomial of 
r(r) with the Lagrange form of the remainder term; that is, we write the curve as 

(4) r(r) = r(0) + r(0)r + ^r(0)r 2 + 1 r (C)r 3 
for some ( £ (0, r) . 

After substituting |4]) into the left-hand side expression in Q and simplifying, 
we obtain that 

(p 4 ,P*+i - ft) = (*± 2 (0) + ^pj/(0) + ipy(0))A^ + Err A 3 , 

where Err is an error term depending on the first three derivatives of r(r), A n and 
i. Note that for any fixed value of i, or more generally if we keep i bounded, this 
error term is bounded. Thus, for any K, if n is sufficiently large and if |z| < K, 
then the sign of this expression is determined by the first three terms. In other 
words, in this case the left-hand side inequality of (J3| is satisfied if and only if 

z± 2 (0) + ip2/(0) + zpy(0) <0 
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Note that k = J^mj- Thus, we have that 

i(pK+l) < ~~PK- 

By a similar argument we obtain that for any value of |i| < K, if n is sufficiently 
large, then the right-hand side inequality of ([3| is equivalent to 

— 1 — -pK < i(pK + 1). 

Thus, for any |i| < K, relint[pi,pi + i] contains an imaginary equilibrium point if 
and only if 

(5) - 1 - ^pn < i(pn+ 1) < -^pn- 

Hence, the number of these edges is no more, and if K is sufficiently large then it 
is equal to, the different values of i that fit in a certain interval of length ^ p ^ +l 

Since the possible values of i are congruent mod 1, this number is either 
or one more. 



\ P K+1\ 



Now we deal with the vertices of P n . Note that pi is an equilibrium point if, and 
only if one of the following holds: 

(6) (-pi,p i+ i -pi) > and(-p i ,p i _i - p^ > 0, 
or 

(7) (-pi,p i+ i - pi) < and(-p i ,p i _i - p^ < 0. 

Using the same kind of argument that we used for edges, we obtain that for suffi- 
ciently large values of n and K and for |i| < K 7 the systems of equation above are 
equivalent to 

1 . . 1 

-pn < i(pn+l) < ~-pn 

and 

1 , . 1 

-^P K < AP K + 1) < 2^ K > 

respectively. 

Clearly, if k > 0, then the first system can be satisfied for no value of i, and 
so is the second one if k < 0. Thus, in both cases, the number of the values of i 



satisfying the corresponding system of equations is either 



pM I 
\ pK +i\] 



i or one more. □ 

If 



4. The equilibria of surfaces 



Similarly like in the previous section, we may assume that the surface, given 
in the form r : D — > E 3 with (u,v) E D = [ui,U2] x [vi,i>2], satisfies the C 3 
property. For simplicity, we examine only the case that the surface has exactly one 
equilibrium point, namely m = r(0, 0) = (0,0, p) with (0,0) € intZ? and p > 0. 
Without loss of generality, we may assume that the tangent plane of r at m is 
horizontal, and that its normal vector n* — 7-„(0,0) x r v (0,0) is an outer normal 
vector of the surface; that is, it points in the direction of the positive half of the 
z-axis. We denote the two principal curvatures of the surface at m by k\ and K2, 
and note that, dealing with generic equilibria, we assume in the rest of the section 
that pni + 1^0^ pn 2 + 1. To be able to define the approximating surface, we 
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assume that r(D) is convex; more specifically, that r(D) is a subset of the boundary 
of a compact, convex set with nonempty interior, containing o in its interior, which 
yields, in particular, that «i,K2 _ 0. 

Consider an equidistant partition of the rectangle D into n 2 homothetic copies 
of ratio -. Similarly like in the previous section, for the indices of the division 
points of the partition we use not necessarily integers, but real numbers that are 
congruent mod 1. More precisely, we set p™ • = r (i U2 ~ Ul , j' " 2 ^" 1 )■ For simplicity, 
for the point p™ •, we may use the notation pij. Furthermore, we set A Ui „ = " 2 ~ Ml , 
A„,„ = ^ and A= 

Let S 1 = {Pij,Pi+i,j,Pi+ij+i,Pij+i}, and set Q — conv(5 U {o}). Since for 
sufficiently large values of n the equilibrium points of P n arc contained in a small 
neighborhood of m and since r u (u, w) and r v (u, v) are continuous functions, we may 
assume that o is a vertex of this polyhedron, and that the rays emanating from o 
and passing through Pi,j,Pi+i,j,Pi+i,j+i,Pi,j+i, respectively, arc in counterclock- 
wise order. 

If S is not coplanar, then either [pij 7 p i+ ij + i] or [p i+ ij 7 pij +1 ] is an edge of Q, 
depending on the position of the points of S. Thus, in this case Q has two triangle 
faces not containing o, namely either \p it j,p i+1 j,p i+li j +1 ] and \pij,pi + i t j + i,pij + i], 
or \pi,j,Pi + i,j,Pi,j+i] and \p i+ltj ,p i+lij+1 ,p itj+1 ]. We denote these faces F^ and 
Ff -. If S is coplanar, we dissect it into two triangles by either \pi,j,Pi+i,j+i] or 
[Pi+i.jiPi,j+i]i an( i regard it as the union of two triangle faces. Now we define the 
approximating polyhedral surface P n as the union of the faces Ff - and Ff - for all 
possible indices The vertices of P n are the points pij, and the edges and the 

faces are those described in this paragraph. We observe that P n is simplicial, and 
that it is not necessarily a subset of the boundary of a convex polyhedron. 

The equilibrium points of P n can be vertices, or relative interior points of edges 
or faces of P n . We assume that each of these equilibria is generic, that is, if q is an 
equilibrium of P n in the relative interior of a face F of dimension 0, 1 or 2, then for 
some neighborhood V of q, the plane passing through and orthogonal to q contains 
no point of V other than those of F. We denote the number of the equilibrium 
points of P n at a vertex p^- satisfying |i| < K and |j| < K, by U n {K). Similarly, 
the numbers of the equilibrium points on an edge or a face of P n , which has a 
vertex p ij:) - satisfying |i| < K and \j\ < K 7 by N n {K) and S n (K), respectively. Our 
aim is to estimate U n (K), N n (K) and S n (K) if n and K are sufficiently large. 

Note that, by the conditions for r described before, we have r„(0, 0) = (x u , y u , 0) 
and r„(0, 0) = (x v ,y v , 0), and that the outer normal vector of the surface at m is 

r„(0,0) x r„(0,0) 



|r u (0,0) x r„(0,0)|| 



= (0,0,1) 



We recall the notions of the first and the second fundamental quantities of a 
surface. In our setting, for these at the point m we have 

E = xl + yl, F = x u x v + y u y v and G = x\ + y 2 v . 

Furthermore, L, M and are the ^-coordinates of r uu (0, 0), r uv (0, 0) and r vv (0, 0), 
respectively. It is well-known in differential geometry that 

(8) EG - F 2 = (r„(0, 0) x r v (0, 0)) 2 = (x u y v - x v y u f, 
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and that 
(9) 



LN - M - , EN - 2FM + GL 

K\K 2 = -frpz FP>") an( i K l + K 2 = 



EG — F 2 ' 1 z EG-F 2 

According to our conditions, we have E, G > 0, -EG — F 2 > 0, L,N < and 
LN - M 2 > 0. 

In our investigation, we examine only the case that AM < L and AM < A 2 iV. 
Since LN — M 2 > 0, these two inequalities are satisfied for some values of A; 
for instance, for the value satisfying L = X 2 N. These conditions ensure that the 
discretization mesh sizes do not differ radically in the u and v directions. We note 
that, in order to obtain the area of a surface by taking the limit of an approximating 
triangulated surface, a similar condition is necessary. 

In the formulation of the main result of this section, we need the following 
notations. 

V 1 =\pL{G + pN), W 1 =\ 2 P N(F + pM), 

V 2 = pL(F + pM), W 2 = XpN(E + pL) 

V 3 = X 2 p(NF — MG), W 3 

V i = \((EG-F 2 )-p(MF-NE)), W t 

V 5 = \(p{NE-MF) + (EG-F 2 )), W 5 



Vq = F + pM, 

V 7 = X 2 p(NF — MG), 

V 8 = X(G + pN), 

V 9 = V 5 , 

V w = Ve, 

Vn - V r , 

v 12 = v s , 

and 



Err,7 



Err s 



4 + 2 



\((EG - F 2 ) - p(MF - LG)), 
p(LF-ME), 

Xp 2 (X\N\ - \M\)(ME — LF), 
W 6 =p{X\N\-\M\)(E + pL), 
W 7 = Xp(\L\ - X\M\){p{LG - MF) + {EG - 
W 8 =p{\L\-X\M\)(F + pM), 
W g = Xp 2 \M\{ME-LF), 
W w = p\M\(E + pL), 
Wu = X 2 \M\(p(LG - MF) ■ 
W 12 = p\\M\(F + pM), 



E max{|K|, \W S \, \V S - W a \, \V S + W s \} 

s=l 

X(EG-F 2 )\( Kl p+l)(K 2 p + l)\ ' 
E max{\V a \, \W S \, \V S - W a \, \V S + W s \} 

s=3 



(EG-F 2 )), 



X(EG- F 2 )\( KlP + l)(K 2 p+l) 



Err 



N 



6 + 



E max{|K|,|W s |,|V; 

s=5 



W S \,\V S + W S \} 



X(EG-F 2 )\(p Kl + l)(pK 2 + l)\ ■ 
Our main result is the following. 

Theorem 2. Let K = (p, K\,K 2 , A) be sufficiently large. Then there is a value n £ 
such that for every n > n e , we have 



U n (K) - 
S n (K) - 
N n (K) + 



P K\K 2 
(KlP+ 1)(k 2 P+ 1) 
1 



(Kip+ 1)(k 2 P+ 1) 

(Kl + K 2 )p 



{Kip+ 1)(k 2 P+ 1) 



< Erry, 

< Err 5 , 

< Errw . 
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We can interpret this theorem in the following way. Our result in Theorem [T] 
can be interpreted in a way that the number of a certain type of equilibrium point 
fluctuates around a specific number that depend only of local quantities of the 
curve at m, within an error one. In general, this fluctuation seems random as we 
compute this numbers for different values of n. Here we have a similar result, with 
the only difference that the error term is not one, but a complicated expression of 
the fundamental quantities of the surface and A. In general, this fluctuation seems 
random; that is, if we compute these numbers for different values of n, then the 
discrepancy between the measured number and the corresponding estimate in the 
theorem 'averages out'. For more information about the sizes of the error terms, the 
reader is referred to Remark [3J and for a more precise treatment of the fluctuation 
to Section [5] The error functions in Theorem [2] provide a global bound and may in 
many cases grossly overestimate the actual errors. In particular, the mesh ratio A 
as well as the topology of flocks may radically influence the actual errors. This can 
also be observed in our example illustrated in Figure [I] showing a triaxial ellipsoid 
with axis ratios a:b:c=l. 23:1. 15:1. 

We start with two lemmas that we need in the proof. 

Lemma 1. Let S C M. d a nonempty Lebesgue measurable set with finite measure 
X(S). Let pk denote the probability that x + S contains exactly k points of the lattice 
1 d , where x is chosen from the cube [0, l] d using the uniform distribution. Define 
the probability distribution a by P(a = k) = Pk- Then the expected value of a is 
E(a) = \(S). 

Proof. First, note that as A(5 I ) is finite, we may write S as the disjoint union of 
finitely many Lebesgue measurable sets with diameters less than one. Let these sets 
be Si, S*2, . . • , S m . Let Oi be the distribution such that P{ai = k) is the probability 
that x + Si contains k lattice points, where x is chosen uniformly from the unit 
cube of M. d . By the linearity of the expectation, we have 

m 

E(a) = Y J E{o- l ). 

i=l 

We show that E(<Ji) = \{Si), from which the assertion readily follows. Without 
loss of generality, we may assume that S C (— 1, 0) d . Then for any x £ [0, l] d , x + Si 
either contains the origin, or it does not contain any point of 7L d . Note that the set 
of points x for which x + Si contains the origin is —Si, and thus, E(<Ji) — A(<Sj). □ 

We remark that, clearly, the lattice Z d in the formulation of Lemma [l] can be 
replaced by any lattice, if x is chosen from the corresponding fundamental lattice 
parallelotope. 

Lemma 2. For every value of K , if n is sufficiently large then for every vertex pi^ 
of P n with \i\ < K \j\ < K we have the following: 

• If M > 0, then E + = [pij,pi + i,j + i] is an edge of P n , and 

• if M < 0, then £L = [p»+i,3 - ,Pi 3+1] is an edge of P n . 

Proof. Recall that if is sufficiently close to m, then o and the endpoints of E + 
and E_ are in convex position, and the rays emanating from o and passing through 
Pi,j,Pi+i.jiPi+i,j+i an d Pi,j+i are in this counterclockwise order in their conic hull. 
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Let S+ denote the sum of the volumes of the tetrahedra [OiPij^i+ijiPi+i j+i] 
and [o,Pi+ij+i,Pi.j+i,Pi } j], and S_ the sum of the volumes of [o,Pi t j,Pi+i.j,Pij+i\ 
and [o,pi+i.j,pi+i J --|-i,pi. :/ +i]. Observe that if S+ — S- > 0, then E + is an edge 
of P n , if 5+ — S- < then E-, and if S+ — S_ = 0, then the two segments are 
coplanar. Note that the signed volume of the tetrahedron spanned by the vectors 
a, b and c is \{a, b X c). Thus, the sign of 5+ — S- coincides with the sign of 

V = (pi,j,Pi+i,j x p i+1 j +1 ) + (Pi + i.j + i,Pi, J+ i x p hj ) 

~(Pi,j,Pi+l,j X Pi,j+l) - (P-i+lj^i+lj + l X Pi,j+l). 

We use the second degree Taylor polynomial of r(u, u) to approximate the surface. 
In other words, we write the surface in the form 

(10) r(u, v) = r(0, 0) + r tt (0, 0)« + r„(0, 0)v + ^r uu (0, 0)u 2 + r uv (0, 0)uv+ 

+ ^r vv (0, 0)v 2 + K^ 3 + K 2 u 2 v + K 3 uv 2 + K 4 v 3 , 

where Ki,K2, K% and K4 are bounded M. 2 —> M 3 functions of u and v. 

After computing V and substituting „ = AA„ „, we obtain a polynomial of 
A u n . In this polynomial, the first nonzero coefficient, which is that of A„ n , is 
My/ EG — F 2 . Thus, for any \i\ < K and \j\ < K, if n is sufficiently large, then 
sign(S' + — S-) = signM. Hence, if M > then E + is an edge of the convex hull, if 
M < then E_, and if M = then it can be both: it is determined by the higher 
degree terms in the approximation. □ 

Now we prove Theorem [2] 

Proof of Theorem [£| We prove the three formulas separately. 

Part 1: The proof of the assertion for U n (K). 

As in the proof of Lemma |2j we approximate the surface r(u,v) by its second 
degree Taylor polynomial. We start with the observation that pi j is an equilibrium 
point if, and only if, each of the eight angles of the form /.(opijPi+SiJ+s ), where 
Si, Sj € { — 1, 0, 1} and (Si, Sj) ^ (0, 0), is acute, or each is obtuse. We note that by 
Lemma [2] it is sufficient to examine six of the angles. 

First, we examine the condition under which /.(opi t jPi+i t j) < \. Note that this 
condition is equivalent to 

( n ) (-Pi,j,Pi+l,j ~Pi,j) > 0. 



After substituting (10 1 and A v n — AA„ „ into this inequality and then expanding, 
the first nonzero term, which does not depend on K\, K2, K3 and K4, is of degree 
two. Since the coefficients of the remaining terms are bounded functions for any 
\i\ < K and \j\ < K, if n is sufficiently large, then (11) is satisfied if and only if 



the first nonzero term is positive. We can rewrite this term using the fundamental 
quantities of the surface, which yields the following inequality: 

(12) X < - l -pL, 
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where X = i(E + pL) + Xj(F + pM). From the condition Z(opijPi-ij) < f, we 
obtain similarly that 

(13) \ P L < X. 

We remark that from Z(opijp i+ ij) > ? and Z(opjj-pj_ij) > | it follows sim- 
ilarly that if n is sufficiently large, then ~\pL < X < \pL. Since L < 0, it is a 
contradiction, and thus, we have that if pij is an equilibrium point, then all the 
eight angles are acute. 

From the condition that the remaining six angles are acute, we have that 

(14) h 2 P N <Y <-h 2 pN, 
where Y = Xi(F + pM) + X 2 j(G + pN), and 



. ±X 2 pN + XpM + ±pL < X + Y < -±X 2 pN - XpM - ±pL 

" \X 2 pN - XpM +\pL < Y-X < -\X 2 pN + XpM - \pL 

Let P denote the convex polygonal region in the (X, F)-plane defined by the 



inequalities in (12), (13), (14) and (15). We determine the shape of P 



Observe first that as L, N < 0, the first four inequalities define a (possibly 



degenerate) rectangle R. Next, we note that the two pairs of inequalities in (15) 



define two (possibly degenerate) infinite strips. Recall that LN — M 2 > 0, which 



implies that the quadratic forms of A on the right-hand sides of ( 15 ) are positive 



semidefinite. Thus, the two infinite strips defined in ( 15 ) are of nonnegative widths 



From (15), it follows also that if M < 0, then the first strip contains R, and 
otherwise the second one. Furthermore, since X\M\ < |L|,A 2 |iV|, the other strip 
contains exactly two (opposite) vertices of R. Hence, P is a (possibly degenerate) 
centrally symmetric hexagon that is obtained as the rectangle R truncated with 
parallel lines at two opposite vertices. This hexagon can be observed in Part (c) of 
Figure [T] 

We have that, for sufficiently large n, the number of the equilibrium points at the 
vertices of P n is equal to the points of a translate of the lattice {(X(i, j), Y(i,j)) : 
i,j € Z 2 } contained in P; or, equivalently, the number of lattice points contained 
in a translate of P. Let Q denote the fundamental parallelogram of the lattice; that 
is, Q = {(X(i,j),Y(i,j)) : i,j e [0,1]}. Then, by Lemmajl] choosing the translate 
of P 'randomly', the expected value of the number of lattice points in it is ^ca(Q) ' 

It is an elementary computation to show that area(P) = X 2 p 2 (LN — M 2 ). On 
the other hand, area(Q) is the absolute value of the determinant D of the linear 
transformation T defining X and Y. Expanding D and using the identities in Q, 
we obtain that 

area(P) p 2 KiK 2 
area(Q) \{pKi + l)(pK 2 + 1)|" 

To prove the assertion it suffices to show that the discrepancy between this ratio 
and the number of the equilibrium points is within Err;/. Because of the properties 
of the expected value, for this we need only show that the difference between the 
numbers of the lattice points contained in two different translates of P is at most 
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Erry, or in other words, that the difference between the numbers of the points of 
Z 2 contained in two different translates of T~ 1 (P) is at most Err^r. 

Note that if we move a convex region parallel to the s-axis, then on each line 
parallel to the rr-axis, the number of the points of Z 2 can change by at most one. 
Thus, the discrepancy is at most the vertical width of the region, rounded up. We 
can observe the same property if we move a convex region parallel to the y-axis. 
Thus, the error term is at most two more than the half-perimeter of the smallest 
axis-parallel rectangle containing T^ 1 (P). Clearly, we can estimate this quantity 
from above by computing the corresponding quantity for the parallelogram T _1 (i?) 
instead of the hexagon T^ 1 (P). Hence, the error term is Erru — 2 + (x max — x m i n ) + 
(2/max - 2/min), where a; max , x min , y max and y min are the extrema of the coordinates 
of the points of T^ 1 (R). Clearly, these extrema are obtained at vertices of T _1 (7?). 
We leave to the reader to verify that by computing these vertices we obtain the 
error term Erru of Theorem [2] 

Part 2: The proof of the assertion for S n (K). 

We prove the statement in three different cases. 

Case 1, M > 0. In this case we need to determine the conditions under which 
there are equilibrium points on [pi,j,Pi+i,j,Pi+i,j+i] or [p,- +1)J - + i,p ij3 - + i,p ij3 -]. 

Consider the face [pij,pi+ij,pi+i.j+i]. First, we note that an outer normal 
vector of this triangle face is 

n = (Pi+l,j-Pi,j) x (Pi+hi+^~Pij) = Pij x Pi+l,j+Pi+l,j x Pi+1 +Pi+l,j+l x Pi,j- 

Our key observation is the following: There is an equilibrium point on the triangle 
if and only if the ray emanating from the origin and perpendicular to the plane of 
the triangle intersects the triangle. In other words, it happens if and only if n is 
contained in the conic hull of the triangle; that is, if the angles between n and the 
inner normal vectors of the faces of the tetrahedron [o,pij,pi+ij,Pi+ij+i], apart 
from —h, are acute. Thus, we have three conditions that determine if there is an 
equilibrium point on the face: 

(n,Pi,j x Pi+ij) > 0; 
(n,Pi+i,j x p l+1 . 3+1 ) > 0; 
(n,Pi + ij+i x p itj ) > 0. 



After substituting (10) and A v ^ n — \A u n into these inequalities and expanding, 
the first nonzero terms, which do not depend on K\,Ki, K3 and K4, are of degree 
4 in A u .„. Expressing them with the fundamental quantities of the surface, we 
obtain that 

(16) XpME<X, \ 2 pMF - \ 2 (EG — F 2 ) < Y, X + Y < XpME + \ 2 pMF, 
where 

X = Xip(LF - ME) + X 2 j(p(MF - NE) - {EG - F 2 )) + ^XpLF - ^X 2 pNE 
and 

Y = X 2 i(p(LG - MF) + (EG - F 2 )) + X 3 jp{MG - NF) + ^X 2 pLG - ^X 3 pNF. 
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Similarly like in Part 1, we estimate the number of the faces with an equilibrium 



point with the ratio of the area of the triangle defined in ( 16 ) to the absolute value 
of the determinant of the affine transformation defining A and Y. An elementary 
computation yields that the determinant D is 

D = X 4 {EG - F 2 ) 2 (nip + l)(n 2 p + 1) 

and, furthermore, that the inequalities in ( |16[ ) define a right triangle T\ the legs of 
which are of length \ 2 (EG — F 2 ). Thus, the ratio of the area of T\ to the area of 
the fundamental parallelogram of the affine transformation is: 

1 



2|(k 1( o + 1)(k 2/ o+1) 



Now we consider the face [pi + i,j + i,pij + i,pi j]. An argument similar to the 
previous one yields the inequalities 
(17) 

X < X 2 (EG- F 2 )- X 2 pMF, Y<-X 3 pMG, ~X 2 pMF - X 3 pMG < X + Y : 

where X and Y are defined as in the previous case. This region is a reflected copy 
T2 of the triangle Ti, and thus for the ratio of its area and the absolute value of 
the determinant we obtain the same quantity. 

Note that as M > and as E + XF + X 2 G is a positive definite quadratic form, 



the open half planes defined with the third inequalities of ( 16 1 and ([17]) overlap. 
Nevertheless, it may happen that T\ and T 2 do not overlap; it happens, for instance, 
when p is sufficiently large. These triangles can be observed in Parts (dO), (dl) and 
(d2) of Figure [IJ 

Case 2, M < 0. We may apply a consideration similar to that in Case 1. The 
faces we need to examine are \pi,j,Pi+i,j,Pi,j+i] and [p l +i, J ,Pi+i,j+i,p l , J +i\- For 
the first face, we have the inequalities 

(18) 0<A, F<0, —X 2 (EG — F 2 ) <Y — X, 
and for the second one, we have 

X < X 2 (EG - F 2 ) + XpME - X 2 pMF, 

(19) X 2 pMF - X 3 pMG - X 2 {EG - F 2 ) < Y, 

Y - X < -XpME + 2X 2 pMF - X 3 pMG - X 2 (EG - F 2 ), 

where X and Y are defined as in Case 1 of Part 2. These two regions define two 
triangles. By computing their area we obtain the same estimate for the number of 
the equilibrium points as in Case 1. 

Case 3, M = 0. In this case the coefficient of the the term of degree four 
in the polynomial V in Lemma [2] does not determine which of the edges E + = 
\Pi,jiPi+i,j+i] an d E- = [pi + i_j,pij + i] belongs to P n . Note that to estimate S n (K) 



we may use the triangles defined by the inequalities in (16) and (17) from Case 1 if 
this edge is E + , and the triangles defined by ( 18 ) and ( |19[ ) if it is £L . Note that in 



both cases, these two triangles degenerate into the same square in the (A, y)-plane 
if M = (and in this case the two triangles do not overlap). This square is defined 

by 

(20) < X < X 2 {EG - F 2 ), -X 2 {EG - F 2 ) < Y < 0. 
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Thus, in the case M = 0, if a triangle with vertices from amongst pij , Pi+ij , Pi+i.j+i 
and Pij+i contains an equilibrium point, then there is exactly one such face of P n 
no matter which edge, E + or E_ belongs to P n . 

We have yet to determine the error term Errs of the theorem for an arbitrary 
value of M . Our method is similar to the one used in Part 1. Let T be the affine 
transformation defined by the formulas for X and Y. Let 5* be an axis parallel 
square of side length X 2 (EG — F 2 ). Then Errs m an y of the cases is at most four 
more than the perimeter of the parallelogram T^ 1 (S). Hence, a simple computation 
yields the required quantity. 

Part 3: The proof of the assertion for N n {K). 

The proof is similar to those on Parts 1 and 2, and thus we just sketch it. 

Case 1, M > 0. We need to find the conditions under which \pi,j,Pi+i,j], 
[Pi,j,Pi,j+i] 01 [Pi,jiPi+i,j+i] contains an equilibrium point. For simplicity we call 
these edges horizontal, vertical and diagonal, respectively. 

Consider a horizontal edge \pij,Pi+i.j]- Note that it contains an equilibrium 
point if and only if the following holds: 

• The angles Z(o,pij,pi+ij) and Z(o,pi + ij,pi,j) are acute. 

• The two angles between the triangle [o,Pij,Pi+ij], and {pij ,pij_i,p i+ ij] 
and [pij,pi + ij,pi + ij + i], are either both acute, or both obtuse. 

Using the second degree Taylor polynomial, for the first nonzero terms we obtain 
the following inequalities: 

(21) 0<X h <E,0<Y h <-h P E(\N + M), 

or 

< X h < E, -XpE(XN + M) <Y h <0, 

where 

X h = -i(E + pL) - Xj(F + pM) - l -pL, 

and 

Y h = Xip(ME - LF) + X 2 j ((EG - F 2 ) + p(NE - MF)) - Xp(XNE - MF). 

Observe that due to our assumption that X\M\ < A 2 |iV|, the second system of 
inequalities is not satisfied for any value of Y h . Note that the first system determines 
a rectangle Rh in the (X, Y) coordinate-system. 

We estimate the number of horizontal edges with an equilibrium point as usual, 
and obtain the quantity 

-Ep(XN + M) 



Z h 



X(EG-F 2 )\(p Kl + l)(p K2 + l)\- 



For the vertical edge \pi,j,Pi,j+i], using similar quantities X v and Y v , we obtain 
another rectangle R v and the estimate 

_ _ -G P (L + XM) 



(EG-F 2 )\(p Kl + l)(p K2 + l) 
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For the edge \pij,pi+ij+i], we have the inequalities 
(22) < X d < E + 2XF + A 2 G, < Y d < XM{E + 2XF + X 2 G) 



0<X d <E + 2XF + AG, XM(E + 2XF + AG) < Y d < 0, 



where 



X d = -i(E + pL + X(F + pM))~j(X(F + pM) + X 2 {G + pN))-^p(L + 2XM + X 2 N), 
and 

Y d = -Xi(p{LF - ME) + X{{EG - F 2 ) + X{LG - MF))) + X 2 j(p(NE - MF)+ 

+(EG - F 2 ) + X(NF - MG)) + ^Xp(2ME - LF + X{NE + 2MF - LG) + X 2 NF). 

As < XM(E + 2XF + X 2 G), we have that only the first system of inequalities, 
determining a rectangle R d , has solutions. Thus, for the number of diagonal edges 
with an equilibrium point, we obtain 

(E + 2XF + X 2 G)pM 
d ~ X(EG-F2)\(p Kl + l)(pK 2 + l)\- 

Thus, for the sum of the numbers of the edges with an equilibrium point on them 
we obtain the estimate: 

-{LN — 2MF + NE)p _ -( Kl + K 2 )p 



(EG~F2)\(p Kl + l)(p K2 + l)\ KpKi + l)(pK2 + l)|' 

The three different types of edges can be observed in Parts (eO), (el), (e2) and 
(e3) of Figure [TJ 

Case 2, M < 0. We need to examine the edges [Pij,Pi+i.j], [Pi,j,Pi.j+i] and 
[Pi+i.jiPi,j+i\- We may apply the approach described in the previous case, and 
obtain similar formulas for the numbers of horizontal, vertical and diagonal edges 
that contain an equilibrium point. In this way, for the sum of the numbers of the 
edges with equilibrium points we obtain the same estimate. 

Case 3, M = 0. Similarly like in Part 2, we may apply the formulas from the 
M > case if E + belongs to P n , and those from the M < case E- does. Note that 
the inequalities that determine if \pij)Pi+i,j] or \p%,j,Pij+i] contain an equilibrium 
point degenerate into the same inequalities for M = 0, no matter which formulas 
we apply. Thus, for Zh and Z v we obtain the same quantities. Furthermore, the 
inequalities that determine if E + contains an equilibrium point in Case 1, and those 
that determine if E— contains such a point in Case 2, determine a planar figure of 
zero area if M = 0, and thus, if M = 0, we have Z d = 0. 

To compute Err at, we may apply the approach used in Parts 1 and 2. We 
remark that if M > 0, then X d = X h + X v and Y d = Y h + Y v , and if M < 0, then 
X d =X h -X v and Y d = Y h -Y v . □ 
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Figure 1. Flocks of an ellipsoid 



Figure [T] shows an ellipsoid with axis ratios a : b : c — 1.25 : 1.15 : 1. 
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• Part aO shows the equilibrium points of the ellipsoid. The stable, 
unstable and saddle points are denoted by si and s 2 , U\ and u 2 , 
and hi and h%, respectively. 

• Part al shows the equilibrium points of P™ near u\ and h±. Faces 
with a stable point are shaded, unstable vertices are marked with 
x , and edges with a saddle point are drawn with bold lines. 

• Part b shows the equilibrium points of P n near h\. The zoomed 
hexagonal region is framed in Part al. 

• Part c shows the unstable equilibrium points near n\ inside the 
hexagonal region P. 

• Part dO shows the stable equilibrium points near hi with the 
triangles T\ and T 2 . These triangles are separately shown in 
Parts dl and d2, respectively. 

• Part eO shows the saddle type equilibrium points near hi with the 
paralelograms R2 and R3. Parts el, e2 and e3 show these 
parallelograms separately. 

The estimates in Theorem [2] yield U n (K) = 20.19, S n (K) = 22.59 and 
N n (K) — 43.78 for the 'average' number of the imaginary equilibrium 
points in the 'vicinity' of hi. For the ellipsoid and the parametriza- 
tion we used in Figure [l] we have U n (K) = 25, S n (K) — 18 and 
N n {K) — 44. The error terms defined in Theorem [2] are the follow- 
ing: Ettu = 15.11, Errs = 30.04, Err^ = 44.98. As we can see, the 
error functions grossly overestimate the actual errors, as we pointed 
out after Theorem [2j 



5. Results on the average of the equilibrium points of P n 

In this section, for x £ R, frac(a;) denotes the fractional part of x. We start with 
the following theorem. 

Theorem 3. We have the following. 

[3Jl Let I C [0, 1] be a closed interval, and let 77 be an irrational real number. 
Let A n = card{frac(fc?7) G / : k = 1, 2, . . . , n}. Then lim ^ exists and is 

71— >CX> n 

equal to the length of I. 
H2 Let X C [0, l] 2 be a closed axis parallel rectangle, and let 771 , 772 be linearly 
independent irrational numbers in the vector space M over Q. Let B n = 
card{( frac(fc?7i), frac(fc?72)) € X : k — 1,2, ...,n)}. Then lim — exists 
and is equal to the area of X. 

In the proof we use the following, well-known result about the simultaneous 
approximation of irrational numbers. A more general statement is proven, for 
example, in Theorem 4.4, on page 45 of [1.6.. 

Lemma 3. Let X be a closed axis parallel rectangle in the unit square [0, l] 2 , and 
let 771 , 772 be linearly independent irrational numbers in the vector space K over Q. 
Then the set {( frac(n?7i), frac(nr72)) € X : n = 1, 2, 3, . . .)} is everywhere dense. 



EQUILIBRIUM POINTS 



17 



Proof of Theorem^ We prove only ([3j2), as the proof of ([3]l) is just its simplified 
version. We set S = [0, l] 2 and p n — ( frac(n?7i), frac(rt77 2 )) . For simplicity, the 
geometric figures in the proof are imagined as factored by S. 

Let a and b denote the lengths of the horizontal and the vertical sides X, respec- 
tively. By the previous lemma, for every S > 0, there is a positive integer N such 
that frac(TV^i) — Si < S and [^c(jVr)i) = $2 < S. Note that the conditions of the 
theorem yield that Si and S% are irrational. 

Let M be the largest value with frac(M7V?72) < 1. Then 



i 



Consider the first MN points. We estimate how many of them are contained in 
X. Pick any value 1 <k<N. We examine the point set P = {Pk,Pk+N>Pk+2Ni •■• iPk+(M-i)N}- 
These points are equidistant points on a line ( mod S) with slope 5% (cf. Figure [2]). 
We count the number of the points of P contained in X. Those that are are located 
on segments such that their projections on the x-axis, possibily apart from the first 
and the last one, are a long. Measured on the x-axis, the distance of the points of 
P is Si, and thus each such segment, possibly apart from the first and the last one, 
contains at least — 1, and at most + 1 points. There are at least ^ — 1 'full' 
segments, and at most + 1 segments of any kind. Hence, the total number of 
points in the rectangle is at least (f^ — lj (j^ — lj and at most (f^ + lj (j^ + ij . 



s 



X 




FIGURE 2 . An illustration for Theorem |3] 

Thus, we have 

^(t- 1 ) fe" 1 ) B M ,N(f x+ l) (A + l) 



< - " i' < 



N 
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Observe that as 61,62 < 6, the inequalities above remain true if we replace 61 and 
62 by 6. 

Clearly, in the way described above we can prove the same inequalities for the 
second MN points, and so on. Thus, for every value of k, we have 

(a-6)(b-6) B kMN (a + 6)(b + 6) 

(1 + 6) 2 kMN (l-<5) 2 ' 

Finally, if kMN < n < (k + 1)MN, then 

B n B (k+1)MN (k + l)MN B (k+1)Mf f < n + NM {a + 6){b + 6) 

n ~ kMN ~ kMN (k + 1)MN ~ n - NM (1 - 5) 2 ' 

and we may obtain similarly that 

n-NM (a-6)(b-6) B n 



n + NM {1 + 5) 2 ~ n 



Now, let e > be arbitrary. Then there is a sufficiently small 5 > such that 
(T~fri^ — + f • Note that MN depends on 6, but it is fixed if 6 is fixed. Thus, 
for the value of 6 > above, there is a positive integer N% such that from n > N% 
it follows that ^5 < ^rrf-, and thus that ^ < ab + e. We obtain similarly 

n+NM ab+ 1 ' n ■' 

the existence of a positive integer N2 such that n > N2 yields that ^ > ab — e. 
Hence, with the notation N* — max{7Vi, ^2}, we have that n > N* implies that 
I — — ab\ < e, or in other words, that lim — — ab. □ 

Remark 1. The interval / and the rectangle X in Theorem [3] can be replaced by 
any Jordan measurable set. 

Theorem [3] and Remark [T] yields the following corollary. 
Corollary 1. We have the following. 

[TJl Let r\ e R \ Q and X C R be a set with Jordan measure X(X). Let 

A n = card(Z n (nr] + X)). Then lim ^ — = X(X). 

[TJ2 Let 771,772 be linearly independent irrational numbers, and let X C E 2 be 
a set with Jordan measure \{X). Let B n = card (Z 2 n (71771,71772) +Xj. 

Then lim ^ — = \(X). 

n— too n 

Proof. Again, we prove only ([l]2). Clearly, we may assume that the diameter of 
X is less than one and that it is contained in [— 1,0] 2 . Then it may contain at 
most one point of the lattice 1? . Observe that (71771,7^772) + X contains a lattice 
point if and only if ( frac(n77i), frac(n?72)) + X contains the origin, which happens 
if and only if ( frac(n?7i), frac(7i772)) is contained in —X. Since the measure of — X 
is X(X), from this the assertion readily follows. □ 

In the previous sections we have seen that the numbers of the different types 
of equilibrium points 'fluctuate' around specific values as we change 71. Roughly 
speaking, in Theorem [4] we prove that if the parameter of the equilibrium point 
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m is positioned 'irregularly' in the parameter range, then the fluctuation 'averages 
out' in the long run. 

Theorem 4. We have the following. 

|3Jl Let r : [ti,ti] — > M 2 satisfy the conditions of Section pj Assume that T ^l Tl 
is irrational. Then, using the notations of Section un if K is sufficiently 
large, we have 



E U k {K) 



lim - — = r— — — 1 — - and 

n-yoo n \pn + 1] 

n 



p\n\ 



lim 

n— > oo 



E s k (K) 

k=l 



1 



n \pn + 1\ 

2J2 Let r : [1*1,7x2] X [ui,«2] I^ 3 satisfy the conditions of Section ^ Assume 
that u "2 Ul an d V2 V - Vl are linearly independent irrational numbers. Then, 
using the notations of Section ^ if K is sufficiently large, we have 



E u k (K) 
,. k=l 

lira 



E N k (K) 



lim — = v ; ' — z/r and 

(->oo n K«iP + l)(ft2P + lj| 

n 

kl 



|(«ip + l)(«2p + l) 



-(«! + K 2 )/0 



lim 

n— > oo 



E S k (K) 

k=l 



1 



|(«ip+l)(/t2p+l)| 



Proof. We prove only (|4]2). Consider, for example, Case 1 of Part 2; that is, the 
case of S n (K) with M > 0. In the proof, we estimated the number of points pij 
such that (i,j) is contained the triangles !T _1 (Ti) and T _1 (T 2 ). Let 771 = ^ "I" 1 and 
^2 = ;f~E^r With the notation used in the proof, we have that frac(i) = frac(nr?i) 
and frac(j) = frac(rwy 2 ). Thus, we need to determine the number of the points of 
the lattice (71771,71772) + in T\ or in T 2 . Hence the assertion follows from the 
previous corollary. In the rest of the cases, we may apply a similar argument. □ 



6. Questions and concluding remarks 

In light of our results in Sections [3] and |4j we introduce the following notions. 

Definition 3. Let m € M. d , where d = 2 or d = 3, be a generic equilibrium point 
of a smooth hypersurface H C M. d with respect to o, and set p = \\m\\. Then the 
imaginary equilibrium indices of m are the quantities: 

• U* = 1 and S* = 1 — ^ttj for d — 2, where k is the (signed) curvature 
|p«+i| |p«+i| v ' 

of H at m; 

. 7T* — IP 2 «l^| AT* _ |/)(ki+« 2 )| 1 c* _ 1 r 

U - \{p Kl +l)(pK 2 + l)\> ^ ~ |( P K 1 + 1)( P K 2 + 1)| D - |(p Kl + l)(p K2+ l)| 101 

d = 3, where n\ and /c 2 are the fundamental curvatures of H at m. 
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In our investigation, we found the equilibrium points of P n with indices in a 
given 'large' neighborhood of or (0,0). It would be convenient to say that if 
we choose a large value of K, then, for a fine discretization, we have found all 
the equilibrium points of P n in this way. Unfortunately, we cannot do that: our 
method does not take into account an infinite sequence p™ k of equilibria if if. is 
not a bounded sequence of Ufe. This happens if, for example, if. is of order ^/rik, 
in which case —> m is still satisfied. This observation leads to the following 
definition. 

Definition 4. Let r be a curve satisfying the conditions in Section|3j with a unique 
equilibrium point to = r(0). If is a sequence of equilibrium points of P n with 

lim if. — oo, then the sequence {p" h } is called an irregular equlibrium sequence. 

fe-voo * 

We may similarly define the notion of an irregular equilibrium sequence for a 
surface. Note that if a curve or surface has no irregular equilibrium sequence, then 
for this curve or surface the estimates for U n (K), N n (K) and S n (K) in Sections [3] 
and [4] hold for the numbers of the different types of all the equilibrium points of 

pn 

A direct computation shows that if r : [r^ra] — > K 2 satisfies the conditions in 
Section [3j and its coordinate functions are polynomials, then r has no irregular 
sequences. This leads to the following question. 

Question 1. Prove or disprove that if the coordinate functions of r are analytic 
functions of r, then the curve has no irregular equilibrium sequence. What about 
surfaces in R 3 ? 

Our method was to approximate our smooth surface S with a polyhedral, and 
thus nonsmooth, surface P n . We found that no matter how close P n is to S, it 
may occur that P n has more equilibrium points than S. This is not the case if we 
approximate S with a surface the first and second partial derivatives of which are 
close to those of S. 

Remark 2. Let r : D — > M? be a C 2 -class surface satisfying the conditions of 
Section |4| Let s > be sufficiently small. Then any g : D — >• W 3 C 2 -class surface 
with the property that the difference between r(u, v) and g(u, v), and between any 
first or second partial derivatives of r and g at (u, v) G D is less, than e, has exactly 
one equilibrium to'. Furthermore, to' is stable, unstable or a saddle point if, and 
only if, to is stable, unstable or a saddle point, respectively. 

Proof. We use the notations from Section |4j 

Let r : D — > M 3 be a C 2 -class surface as in Section |U Note that this surface 
has a stable, unstable or saddle point at to = r(0, 0) if exactly two, zero and one 
of pK\ + 1 and pK2 + 1 is positive. Thus, since r is C 2 , if we approximate r by 
g in a way that the surface, and its first and second partial derivatives change at 
most e for some small e > 0, then the signs of pn\ + 1 and pK-2 + 1 do not change 
in a neighborhood of (/(0, 0). Clearly, if e is sufficiently small, g has no equilibrium 
outside this neighborhood, and the type of any equilibrium inside is the same as 
that of to. Hence, our remark follows from the Poincare-Hopf theorem. □ 
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Remark 3. In Section |4j to prove Theorem [2] we made the assumptions that 
A|M| < \L\ and X\M\ < X 2 \N\. Nevertheless, in the opposite case, a slight modi- 
fication of the proof of the theorem shows that the estimates in the theorem are, 
apart from the error terms, lower bounds for U n (K), N n (K) and S n {K). 

Remark 4. We have not examined how large the error terms in Theorem [2] are 
compared to the estimates of U n (K), N n (K) and S n (K). Note that the estimates 
are approximated by the areas of certain convex regions, and that the error terms 
are approximated by the sum of the half perimeters of certain parallelograms cir- 
cumscribing these regions. Thus, roughly speaking, if the estimates are 'large', then 
'probably' the error terms are 'relatively small'. 

Remark 5. In Theorem |4j we described the number of equilibrium points in a 
more precise way than in Theorems [T] and [2] under a special condition determined 
by the location of the parameters of the equilibrium point of r in the domain of r. 
Note that choosing the parameter range 'randomly', the conditions of Theorem [4] 
are satisfied with probability one. 

Remark 6. Although we deal only with generic equilibrium points, Theorem [2] 
also gives a prediction for the degenerated cases as up —> —1. Assuming p=l and 
Ki < K2, as Ki — > — oo and k-i — > — 1 we obtain the following limits: 

(23) lim lim U* = lim lim U* = oo 

Ki — ► — OO K-2 — >—l > — 1 Ki — — OO 

(24) lim lim N* = lim lim N* = oo, 

Kl — y — OO K2 — ► — 1 K2 — Hi Kl — ^ — OO 

however, the limit of S* is undefined: 

(25) lim lim S* = 0, 

K,2 — > — 1k± — > — OO 

(26) lim lim S* = oo. 

Kl— > — OO K2— 1 

Degenerated flocks corresponding to this limit are illustrated on a flattened oblate 
spheroid in Parts (a) and (b) of Figure [3] In accordance with (23)- (24) we can 



observe large numbers of unstable and saddle points, however, the number of stable 
equilibrium points may be low (cf. Equation ( [25] ) and Part (a) of Figure [3]) or high 
(cf. Equation (26) and Part (b) of Figure [3|. Similarly, as K\ — » — 1 and k 2 — > 0, 



S* and N* converges to oo, however, the limit of U* is undefined. This situation 
is illustrated on an elongated prolate spheroid in Parts (c) and (d) of Figure [3j 

Remark 7. Imaginary equilibrium indices were computed on natural pebble sur- 
faces and were compared with the real number of equilibrium points (Figure [4]). 
Pebbles were digitized by a high-accuracy 3D scanning method, resulting a dense 
point cloud. Since equilibrium points are located on the convex hull of the surface, 
we identified equilibria on the convex hull of the point cloud [10] . Imaginary equi- 
librium indices were computed by using a curvature estimation method [S] on the 
triangular mesh. Although the convex hull of a pebble is clearly not an equidis- 
tant discretization, as Figure [4] illustrates, imaginary equilibrium indices give a fair 
estimate on the actual (integer) number of equilibrium points. 

Clearly, imaginary equilibrium indices can be defined for hypersurfaces H C K d 
for d > 3 as well. Hence, a natural question is to ask about higher dimensional 
analogues of our theorems. 
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Figure 3. Flocks on an oblate and on a prolate spheroid 




Figure 4. Flocks on a pebble. Imaginary equilibrium indices 
(S* ,U* , N*) predicted by ^ and actual numbers of equilibrium 
points (S, U, N) measured on the pebble: 
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